Observation of Two-source Interference in the Photoproduction Reaction 
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In ultra-peripheral relativistic heavy-ion collisions, a photon from the electromagnetic field of one 
nucleus can fluctuate to a quark-antiquark pair and scatter from the other nucleus, emerging as 
a p°. The p° production occurs in two well-separated (median impact parameters of 20 and 40 
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fermi for the cases considered here) nuclei, so the system forms a 2-source interferometer. At low 
transverse momenta, the two amplitudes interfere destructively, suppressing p° production. Since 
the p° decays before the production amplitudes from the two sources can overlap, the two-pion 
system can only be described with an entangled non-local wave function, and is thus an example 
of the Einstein-Podolsky-Rosen paradox. We observe this suppression in 200 GeV per nucleon-pair 
gold-gold collisions. The interference is 87% ± 5%(stat.) ± 8% (syst.) of the expected level. This 
translates into a limit on decoherence due to wave function collapse or other factors, of 23% at the 
90% confidence level. 



Relativistic heavy ions carry strong electromagnetic 
fields which can be treated as sources of quasi-real virtual 
photons. When two ions collide, a large variety of two- 
photon and photonuclear interactions can occur [l|. In 
coherent vector meson photoproduction, a photon from 
the field of one nucleus fluctuates into a virtual quark- 
antiquark pair which scatters elastically from the other 
nucleus, emerging as a real vector meson. p° photopro- 
duction has a large cross section, 8-10% of the hadronic 
cross section for gold-gold collisions at a center-of-mass 
energy of 200 GeV per nucleon-pair 0, H, 3 . Photopro- 
duction can occur at large impact parameters, b. For p° 
photoproduction the median b is about 46 fm 

The (qq)N scattering that produces p° occurs via the 
short-ranged strong force; the p° is produced within one 
of the two ions. The p° source consists of two well- 
separated nuclei. There are two possibilities: either nu- 
cleus 1 emits a photon which scatters off nucleus 2, or vice 
versa. These two possibilities are indistinguishable, and 
are related by a parity transformation. Vector mesons 
have negative parity, so the two amplitudes combine with 
opposite signs. The nuclear separation can be accounted 
for with a transverse momentum (pr) dependent phase 
factor. The cross section is [61 



A(j>t, b, y) - A(p T , 6, -y) exp (ip T ■ b) 



(1) 

where A{px, b, y) and A(j>t, b, —y) arc the amplitudes at 
rapidity y for p° production from the two photon direc- 
tions. We take % = c = 1 here. At mid-rapidity the 
amplitudes for the two directions are equal, and 



ct(pt,6,0) = 2|A(pt,6,0)| ; 



1 — cos (pr ■ b) 



(2) 



The system acts as a 2-slit interferometer, with slit sepa- 
ration b = \b\. The cross-sections at different b are added, 
and the px spectrum is obtained by integrating Eq. (1) 
over b. p° production is suppressed for pr ^ 1/(6), where 
(b) is the mean impact parameter. 

The p° rapidity y and mass my and the photon en- 
ergy ki are related by &i,2 = (my/2) exp(±y) where the 
subscript refers to the two directions. Away from y = 0, 
ki ^ k%, so A(pT,b,y) ^ A(pr,b, —y), and the interfer- 
ence in Eq. ([I]) is less than maximal. 

There are two theoretical calculations of this interfer- 
ence. Klein and Nystrand [f| calculated the interference 



using a detailed nuclear form factor, averaging the pho- 
ton flux over the nucleus. Hencken, Baur and Trautmann 
used a more detailed model of the photon profile and 
a Gaussian form factor for the nucleus Q. This work 
only considered production at mid-rapidity (y = 0), and 
so cannot be directly compared with the data presented 
here. At y = 0, the two calculations agree quite well. 

If the p° production phase depends on the photon en- 
ergy, this would introduce a y— dependent phase shift into 
Eq. fl}. This is not expected in the soft-Pomeron model 
[8j, and we assume that this phase difference is negligible. 

The produced p°s decay almost immediately at two 
well-separated points, so any interference must develop 
after the decay, and involve the 7r + 7r~ final state. Since 
the pions go in different directions, this requires an en- 
tangled 7r + 7r _ wave function which cannot be factorized 
into separate 7r + and n~ wave functions; this is an ex- 
ample of the Einstein-Podolsky-Rosen paradox [tj [l(| • A 
measurement of the two-source interference is sensitive to 
any loss of quantum mechanical coherence, be it due to 
interactions with the environment l_l| or as a character- 
istic of the p° decay. 

Interference is also expected when the p° photoproduc- 
tion is accompanied by mutual Coulomb excitation of the 
two nuclei. This reaction proceeds primarily via three in- 
dependent single-photon subreactions (one to excite each 
nucleus, and one to produce the p°) At a given 6, the 
cross-section for the subreactions factorizes; the probabil- 
ity for an n-photon reaction is P n {b) = EliLi Pi{b) 0' 
where Pj (b) is the probability for subreaction i. There- 
fore, these multi-photon reactions have much smaller 
(&) and the effect of interference extends to higher px 
[5(. Because of the different (b), multi-photon interac- 
tions are important for studying this interference. The 
Klein-Nystrand model uses measured photonuclear cross- 
sections for the mutual Coulomb excitation [l3T |. while 
Hencken, Baur and Trautmann used the Giant Dipole 
Resonance, plus a correction. For e + e~ production ac- 
companied by nuclear breakup, using the measured mu- 
tual breakup cross-sections rather than the Hencken Baur 
and Trautmann approach leads to a 20% larger cross- 
section 14 1; a similar difference may apply for p° photo- 



production. 

In this letter we measure two-source interference in 
200 GeV per nucleon-pair gold-gold collisions by study- 
ing the transverse momentum {pr) spectrum of photo- 
produced p°s. These data were taken with the STAR 
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detector. The major detector component used here is a 
central Time Projection Chamber (TPC) [15] in a 0.5 T 
solenoidal magnet. The TPC tracked charged particles 
with pseudorapidity \rj\ < 1.0. We used two trigger detec- 
tor systems, the Central Trigger Barrel (CTB) and two 
Zero Degree Calorimeters (ZDCs). The CTB consisted 
of 240 scintillator slats surrounding the TPC, detecting 
charged particles with pseudorapidity \rj\ < 1.0 The 
ZDCs detected neutrons emitted by the dissociating gold 
nuclei with virtually unc hang ed longitudinal momentum 
(100 GeV/c) and low p T (nf. 

Data were collected with two different trigger condi- 
tions. The first was a topology trigger which selected 
events with roughly back-to-back pions in the CTB 
Nearly vertical pairs were excluded, to reduce contamina- 
tion from cosmic rays. The second, minimum bias (MB) 
trigger selected p° accompanied by mutual dissociation. 
In these events, both nuclei broke up and released neu- 
trons into the two ZDCs. The cross-section for p° pro- 
duction accompanied by mutual Coulomb excitation is 
about 7% [j| of the total p photoproduction cross-section, 
so the two datasets are essentially independent. 

Events were required to have net charge zero and ex- 
actly two reconstructed tracks which formed a vertex less 
than 50 cm longitudinally from the center of the TPC for 
the MB sample, and 100 cm for the topology sample. The 
difference is because events from the CTB based trig- 
ger were distributed more broadly along this axis. For 
the topology data, we exclude events with \y\ < 0.05 to 
remove the remaining contamination from cosmic rays, 
where a single muon track could be reconstructed as two 
tracks with net charge 0, pr = and y — 0. All tracks 
were assumed to be pions, and were required to have a 
7T7T invariant mass 550 MeV/c 2 < M W7r < 920 MeV/c 2 . 
These criteria produced a clean set of p° events, at some 
cost in efficiency. The 550 MeV/c 2 mass cut removes 
background from misidentified two-photon production of 
lcpton pairs. The background, estimated from the like- 
sign pion pairs, was small, 1.4%. 

To understand the effect of detector resolution, p° 
events were generated following the Klcin-Nystrand dis- 
tributions, and passed through the detector simulation 
and reconstruction. Figure [1] compares the rapidity and 
distributions of the data and simulations. The 
agreement is good for the minimum bias data, less so 
for the topology data. This is most likely due to an im- 
perfect topology trigger simulation; the effect of this is 
treated as a systematic error. 

The applied cuts select both directly produced ir + n~ 
pairs [18| and p°. Direct n + n~ and p° decays are indis- 
tiguishable, so the two processes interfere. The p° mass 
peak and direct pion fraction are consistent with ear- 
lier gold-gold photoproduction studies . The two sub- 
channels should have the same spin/parity and quantum 
mechanical behavior, so we do not distinguish between 
them. If tt + tt~ pairs are produced by a different produc- 
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FIG. 1: Rapidity (left) and M„ (right) of the tt+tt" distribu- 
tions for the topology (exclusive p°, top) and MB (Coulomb 
breakup, bottom) samples. The points with statistical error 
bars are the data, and the histograms are the simulations. 
The 'notch' in the topology data around y = is due to the 
explicit rapidity cut to remove cosmic-ray backgrounds. 



tion mechanism, such as Pomcron-Odderon interactions 
[lij . then this assumption might not hold. 

Figure [5] compares the uncorrected < 0.5 MB data 
and simulations based on Refs. with and without 

interference, as a function of t± = p\. At RHIC energies, 
the longitudinal component of the 4-momentum transfer 
is small, so t± « t. The measured dN/dt spectrum is 
roughly exponential, but with a significant downturn for 
tj_ < 0.0015 GeV 2 , consistent with the predicted inter- 
ference. The no-interference histogram is almost expo- 
nential, dN/dt oc exp (—kt±), where k is related to the 
nuclear radius 0, |2Q|, even though the Klein- Nystrand 
calculation uses a Woods-Saxon distribution for the gold 
density. The Hencken-Baur-Trautmann calculation uses 
a Gaussian distribution for the nuclear density, but is 
also fairly well fit by an exponential. The interference in 
different y ranges is determined using a Monte Carlo sim- 
ulation which follows the Klein-Nystrand calculations. 

Figure [3] shows the efficiency corrected MB and topol- 
ogy data. All four panels show a dip as t± — > 0. As 
expected, this dip is broader for the MB data because 
(6) is smaller. The suppression at t± = is larger for 
the small-rapidity samples because the amplitudes for the 
two photon directions are more similar. The efficiency is 
almost independent of px, so Fig. [2] is not very different 
from the efficiency corrected t± spectra in Fig. [3] The 
main effect of the detector response is px smearing due 
to the finite momentum resolution. 

The dN/dt spectrum is fit by the 3-parameter form 

dN 

^=Aexp(-kt)[l + c(R(t)-l)], (3) 



where 
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FIG. 2: Raw (uncorrected) p° t^-spectrum in the range 
0.0 < \y\ < 0.5 for the MB data. The points are data, with 
statistical errors. The dashed (filled) histogram is a simula- 
tion with an interference term ("Int"), while the solid his- 
togram is a simulation without interference ("Nolnt"). The 
handful of events histogrammed at the bottom of the plot are 
the wrong-sign (tt + tv + + tt~-k~) events, used to estimate the 
combinatorial background. 
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FIG. 3: Efficiency corrected t± spectrum for p° from (top) 
minium bias and (bottom) topology data, for mid-rapidity 
(left) and larger rapidity (right) samples. The points are the 
data, while the solid lines are the results of fits to Eq. (J3j . 



is the ratio of the simulated tx _s P 0C t ra with and without 
interference. For t± > O.OlGeV 2 , R(t ± ) — > 1, but for 
i_L < 0.01 GeV 2 , R(t ± ) ^ 1. A is the overall (arbitrary) 
normalization and c gives the degree of spectral modifica- 
tion; c = corresponds to no interference, while c = 1 is 
the predicted Klein-Nystrand interference. Table 1 gives 
the fit results. 

R(t±) was determined using a simulation that in- 
cludes the detector response, and then fit to two ana- 
lytic functions: R(t±) = T,f =0 ai/(t± + O.O^GeV 2 ) 4 and 
R(t±) = £f =0 aiix- Our results use the first polynomial 
with n = 5; the second polynomial and different values 
of n were used to estimate the fitting uncertainties. 

The weighted average of the four c values is c = 
0.84 ± 0.05. The k values for the MB and topology 
datasets differ by 15%. This may be due to the differ- 



Dataset A k c xV 

(GeV 2 ) DOF 

MB, \y\ < 0.5 6, 471±301 299 ± 12 0.92 ± 0.07 45/47 
MB, 0.5< \y\ < 1.0 5,605±330 303 ± 15 0.92 ± 0.09 76/47 
T, 0.05 < \y\ < 0.5 11,070±311 350 ± 8 0.73 ± 0.10 53/47 
T, 0.5 < \y\ < 1.0 12, 060±471 333 ± 11 0.77 ± 0.18 64/47 



TABLE I: The results of fitting Eq. |3]to the four data sets. 
Here, T is for topology. The /DOF are discussed in the 
text. 

cnt b distributions. The photon flux scales as 1/b 2 , so 
the photon flux on the 'near' side of the nucleus is larger 
than on the 'far' side. As b decreases, p° production is 
increasingly concentrated on the near side, and the ap- 
parent production volume drops, reducing k. A calcula- 
tion with different assumptions may predict a different 
electric field variation which leads to a smaller difference 
in k. 

Two of the fits have x 1 /DOF significantly larger than 
1 . The \ 2 did not decrease with different fit functions for 
R(t), variations of the nuclear radius in the interference 
calculations, background level, or modifications to the 
detector simulation. When the /DOF > 1, we scale 
up the fit errors on c by yjx 2 ' I 'DOF ; this excess error 
may have theoretical and/or experimental origin. With 
the scaled errors, the weighted average is c — 0.86 ±0.05. 

Systematic errors come from instrumental effects, 
background, fitting, and theoretical issues. The major 
instrumental effects were due to the topology trigger; we 
apply a 10% systematic error to the topology data to 
account for this. 

This analysis is sensitive to any p° py-dependent effi- 
ciency variation. The decay pions have a typical px of 
about 300 MeV/c, where the detection efficiency is high 
and almost p^-independent [TH ]. However, the p° pt res- 
olution, about 7.5 MeV/c, smears the t± spectrum in 
the two lowest t± bins. To study detector effects, we fit 
the raw (uncorrected) t± spectrum with the raw Monte 
Carlo output; this reduced c by 18% [2l|, mostly due to 
the pt smearing. We assume conservatively that the de- 
tector simulation is only 80% effective, and assign a 4% 
systematic error on c to account for non-trigger detector 
effects. 

Backgrounds were estimated by including like-sign 
pairs (7r + 7r + + 7r~7r~) in the fits, c changed by less than 
0.5%. We assign a 1% systematic error due to back- 
grounds. 

The uncertainty due to fitting was evaluated by com- 
paring results using the two different polynomial forms 
of R(t) for both n = 4 and n = 5; c varied by an average 
of 1%. The effect of an imperfect form factor model was 
studied by varying the nuclear radius in the simulations. 
A ± 20% change in nuclear radius changed c by 3%. We 
assign a 4% systematic error due to the fitting procedure. 
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The theoretical uncertainties are difficult to evaluate. 
Our simulation follows Refs. [B|, @| in detail, but those 
calculations themselves contain uncertainties. The two 
theoretical models agree well for exclusive p produc- 
tion. For p° production accompanied by mutual Coulomb 
excitation, there is some disagreement, but the Klein- 
Nystrand model has a more detailed excitation calcula- 
tion, and so may be more accurate. We assign a 5% 
systematic error due to theoretical issues. 

Combining these systematic errors in quadrature re- 
sults in an 8% (13%) systematic error for the MB (topol- 
ogy) data. Adding the four results in quadrature, includ- 
ing the systematic errors, leads to an interference that is 
87 ± 5(stat.) ± 8%(syst.) of the expected level. 

Because p°s decay so rapidly, 7/3ct <C (b) , the p° decay 
points are well separated in space-time, and the two am- 
plitudes cannot overlap and interfere until after the decay 
occurs. The interference must involve the tttt final states 
0. This interference is only possible if the post-decay 
tttt wave functions retain amplitudes for all possible p° 
decays, at least until the wave functions from the two ion 
sources overlap. The tt + tt~ wave function is not factor- 
izable and is thus an example of the Einstein-Podolsky- 
Rosen paradox [l(| ■ Unlike previous tests of non- locality, 
the interference involves continuous variables, momen- 
tum and position 0]. 

In conclusion, we have measured the interference be- 
tween p° production at two sources (the two nuclei) by 
observing the tt + it~ decay products. We observe the in- 
terference at 87 ± 5(stat.)±8 (syst.)% of the expected 
level. This shows that the final state wave function re- 
tains amplitudes for all possible decays, long after the 
decay occurs. The maximum decoherence (loss of inter- 
ference) is less than 23% at the 90% confidence level. 
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